Probing the pseudogap phase of cuprates by a Giaever transformer 
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We develop a theory of the rectification effect in a double-layer system where both layers are 
superconductors, or one of the layers is a normal metal. The Coulomb interaction is assumed 
to provide the dominant coupling between the layers. We find that superconducting fluctuations 
strongly enhance the drag conductivity, with rectification most pronounced when both layers are 
superconductors. In view of their distinct dependence on temperature near T c and layer separation, 
drag measurements based on a Giaever transformer could distinguish whether rectification occurs 
due to fluctuating pairs or inductively coupled fluctuating vortices. 
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Introduction- The origin of the pseudogap phase in 
high temperature cuprate superconductors is the subject 
of on-going debate [1] . One possible explanation of this 
phase is the presence of preformed Cooper pairs which de- 
plete the density of states near the Fermi level [2] • This is 
consistent with photoemission studies which reveal that 
the Fermi surface breaks up into disconnected arcs [3] 
that should have profound implications for the nature of 
transport in the pseudogap phase [4]. 

In this connection, recent observations [6] of a large 
Nernst signal in the pseudogap phase has generated much 
interest. Two competing theories may account for this 
phenomenon. On the one hand, in proximity to the su- 
perconducting transition temperature T c , it is natural to 
attribute a large Nernst response to fluctuating vortex- 
like excitations which carry entropy. On the other hand, 
there exists a close correspondence of the Nernst signal 
with the nuctuational diagmagnetism [7 , which points 
towards the Nernst effect as originating from fluctuat- 
ing pairs. Existing theories based on Ginzburg-Landau 
or diagrammatic approaches [8HIU] give a good descrip- 
tion of the Nernst data. It is thus desirable to have an 
additional probe of the pseudogap phase which can dis- 
tinguish transport contributions arising from fluctuating 
vortices or preformed pairs. In this paper, we suggest 
using a Giaever transformer |11) for such a probe and 
develop the theory for this. 

Measurements on double-layer systems [12] have pro- 
vided us invaluable information about the scattering 
mechanisms underlying transport. A typical experiment 
consists of driving one of the subsystems (active layer) 
and measuring an induced response in the other one (pas- 
sive layer). Physically, by driving a current in the ac- 
tive layer, one creates an electromagnetic environment to 
which the passive layer is exposed. These fluctuations are 
rectified, thus producing a finite drag response. Under 
the condition of no tunneling between the layers, momen- 
tum can be transferred only via electron-electron scatter- 
ing due to the inter-layer Coulomb interaction. In con- 
trast, if vortices are present, then inductive coupling be- 
tween the layers becomes important. Experimentally, the 



measured quantity is the drag resistance, Rd — Vd/I, 
which is the ratio between the induced voltage Vd in the 
passive layer to the driving current I in the active one. 
The temperature dependence of Rd (T) is extremely sen- 
sitive to the microscopic properties of the system. For 
example, in Fermi liquids Rd cx T 2 , which results from 
the phase space available for scattering. In contrast, in 
disordered systems Rd oc T 2 lnT [1 3j . while for non- 
Fermi liquids, like composite fermions in quantum Hall 
bi-layers, the drag exhibits a very different temperature 
dependence, namely Rd oc T 4 / 3 [11]. Many other inter- 
esting examples can be found in the literature |15j . 

In recent work [T^] , the rectification effect was studied 
in the context of a magnetic field driven superconductor- 
insulator transition and attributed to fluctuating vor- 
tices. With some modifications, this theory can be ap- 
plied to the case of cuprates under consideration here. 
However, we take the point of view that the dominant 
contribution to the rectification signal in the pseudogap 
phase is due to fluctuating pairs. Below, we develop the 
corresponding theory within a linear response Kubo for- 
malism. 

Formalism- We consider both a symmetric supercon- 
ductor bi-layer and a non-symmetric one where one layer 
is replaced by a normal metal. The building blocks of 
the theory are the normal and superconductor Greens 
functions 
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In what follows we assume a free electron dispersion rela- 
tion, £ p = p 2 /2m*— £p, which in two dimensions implies a 
constant density of states, simplifying the energy integra- 
tions. We assume a momentum dependent d-wave pseu- 
dogap of the form A p = [cos(p x ) — cos(py)] ~ Acos2$ p . 
In Eq. |l]) £„ (s) = £ + 7„( s ) sign(e), with ~f s ~ T being the 
scattering rate observed by photoemission in cuprates, 
and "f n = l/2r„ a constant scattering rate assumed for 
the normal metal. 

For the fluctuating Cooper pairs in a superconductor, 
we take the standard expression for the propagator whose 
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retarded component has the form |17j 

L R (Q,n) = - 1 1 
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n = 23 accounts for electron-hole asymme- 
try, N s is the density of states, D s the diffusion constant 
and e = ln(T/T c ). 

We assume a screened Coulomb potential as the lead- 
ing source of coupling between the layers 
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Here we have neglected retardation effects, d is the 
spacing between the layers and n n r s \ — 47re 2 7V n ( s ) is 
the Thomas-Fermi screening momentum in the normal 
metal (superconductor). Applying a linear response anal- 
ysis [13] to a Giaever device (Fig. [iji), we have for the 
drag conductivity 
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FIG. 1: (Color online) (a) Schematic representation of the 
rectification effect due to preformed Cooper pairs in a double- 
layer superconductor Giaever device. Diagrams (b) and (c) 
depict the rectification coefficients r n (s) i n the normal metal 
and the superconductor, respectively. 



where r ra ( s ) is the nonlinear susceptibility (rectification 
coefficient) within each layer. Knowing <j£> and invert- 
ing the conductivity tensor, one finds the drag resis- 
tivity pd = &D/{&nVs), and finally the drag resistance 
Rd = Pd((-/w), with £ and w being the length and width 
of the two-dimensional strip. In the normal metal, the 
rectification coefficient is given by (Fig. [T|d) 

T n (q,iu)i,iuj2) = v p G n (p±,iei) 

P,ei,± 

G„(p T ,iei±iuj 1 )G n (p T ,iei±iuj2) (5) 

with p± = p ± q/2, and e t = 2ttT(1 + 1/2) being 
the fermionic Matsubara frequency. The susceptibility 
T n (q,u>) which enters Eq. Q is obtained from Eq. ^ 
upon appropriate analytic continuation icJi,2 — > u; ± iO 
such that T n (q,u) = T RA (q,uj,uj), where indices R(A) 
stand for the retarded (advanced) components in the 
complex frequency plane. In the superconducting layer 
we have (Fig. [lj) 

T s (q,iwi,iu) 2 ) = 4eT B(Q±,iu lt2 , ift m ) 

A(Q,q,iQ m + iui,iSl m )A(Q,q,iQ m ,ito m + ioj 2 ) 
L(Q T ,iQ m )L(Q±,iCl m ± iu>i)L(Q± , iQ m ±«w 2 ) (6) 

with Q± = Q ± q/2 and Q m = 2^mT being the bosonic 
Matsubara frequency. The scalar (Coulomb) vertex func- 
tion defining T s has the explicit structure 

A(Q,q,itt m +iuii,iQ m ) = T^G s (p,iei) 

P,Sl 

G S (Q -p,iCl m + iwi - isi)G s (p - q,ie t - iujx) (7) 



while the current vertex reads 
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Similar to the previous case, T s (q,ui) entering Eq. Q 
is obtained from Eq. Q after analytic continuation to 
the real energy axis iu)x,2 — >■ w ± iO. From Eq. (|6| one 
sees an analogy between rectification and the Nernst ef- 
fect. Indeed, the Nernst is given by a Feynman graph 
of three currents connected by fluctuation propagators 
similar to Fig. [IJ;. The rectification coefficient f s has 
the same structure with the only difference being that 
two of the current vertices are replaced by a scalar field 
which couples to the electromagnetic potential in the 
other layer. This point is illustrated in Fig. [T^,, which 
schematically represents the rectification effect due to 
fluctuating pairs. In order to find ctd, we need an ex- 
plicit form for r„( s ) (q, u>) which we derive below. 

Results- We start from Eq. ^ and convert the Mat- 
subara sum over the fermionic frequency into the contour 
integral TJ2 Sl -> § c Mi tanh ^G n (e)G n (e±tuj 1 )G n (e± 
1LO2), where the momentum variables of the Greens func- 
tions have been suppressed for brevity. The contour is a 
circle with three branch cuts at Im(e) = {0, — U}%, — Wa}- 
The outer branches of the contour that have all three 
Greens functions of the same causality do not contribute. 
The remaining parts, after analytic continuation, give us 
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in the ballistic case 
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where we assumed f^°° de [ tanh — tanh^] = 
2w (low temperature limit). We next expand the 
Greens function to linear order in q and perform 
the remaining momentum integration by approximating 
E P #p(#p -£)••• = ^p^^ ■ • • which accounts 

for the curvature of the electronic dispersion. Since 
max{w,g 2 /2m*} <C £f we find to leading order 
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which agrees with earlier calculations [13] ■ One impor- 
tant remark here is that for the linearized dispersion, 
= v fP, rectification vanishes due to electron- hole sym- 
metry. 

We proceed with the derivation of T g (q,u}), taking 
into account a temperature dependent scattering rate 
7 S ~ T appropriate for cuprates [Eq. Q]. Starting 
from Eq. ([6| we convert the bosonic Matsubara sum into 
the contour integral T J2n m ~* j^i £ dficoth + 
iui\, Q)A(0, £l+iui 2 ) L(fl) L((l±iuJi) L(tt±iuj2) with three 
branch cuts at Im(f2) = {0, T^i, ^^2}- After analytic 
continuation and expansion of the fluctuational prop- 
agator [Eq. pj)] to leading non-vanishing order in the 
electron-hole asymmetry Tq, we find 
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When deriving 
this expression, we made use of the following key ob- 
servations, (i) Near T c the characteristic bosonic ener- 
gies are smaller than the fermionic ones {DQ 2 , il, u} ~ 
T-T c <r. This allows us to compute the vertex func- 
tions [Eqs. ([7|-([8])] by neglecting their dependence on the 
bosonic modes. This step is legitimate since both A and 
B contain a fermionic loop with e ~ T . (ii) Under the 
above approximation, one then finds from Eqs. ([7| and 
^ with the help of Eq. Q 
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where scaling functions 
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account for the suppression of the vertex due to the pseu- 



dogap [ID]. In Eq. (14), £ = Hv F /nA and i] = ttD/8A. 



(iii) We also exploited the fact that Q ^ q and set 5 — >• 
in the propagators since rectification is determined by 
small momentum transfer between the layers, namely 
Dq 2 ~ (T-T c )/(K n(s) d) < T-T c ~ DQ 2 . (iv) Our final 
comment regarding Eq. (11) concerns the electron- hole 



asymmetry. In general Tn enters both the fluctuational 
propagator L and the scalar vertex A which, in fact, is 
a consequence of gauge invariance [13l [18] . We find that 
the most singular in T~T C contribution to up comes from 
keeping Tq in L while taking the scalar vertex as the bare 
one. The technical reasoning is that by expanding L in 
To, one effectively raises the power of the fluctuational 
propagator that in the end translates into an extra factor 
of e in an after energy integrations. Conversely, keeping 
Tq in A would also lead to a nonvanishing contribution 
to T s , but with a subleading dependence on T — T c near 
the transition. 



We can now simplify the integrand in Eq. (11) even 



further by exploiting the separation of the energy scales 
{O, w} <C T, thus approximating J"(fi,o;) ~ 11 ~ 
and using Eq. ^ which gives us 
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Here we introduced dimensionless variables = ttJI /8T 
and Aq = ttD s Q 2 /8T+e. The remaining energy and mo- 
mentum integrations can be completed analytically with 
the final result 
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where the Ginzburg-Landau time tql — ir/(8Te) was in- 
troduced. By combining Eqs. (10) and (16) in Eq. Q, 
we find the drag conductivity due to fluctuational rectifi- 
cation by preformed Cooper pairs in the non-symmetric 
normal metal-superconductor transformer (restoring K) 
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where a Q = £, \ ns (d) = { ^% , f(T) = (£) V/3 
. In the immediate vicinity of the transition, 



and r, = 



T — T c <C T c , where f(T) is constant, the rectification 

has a strong power-law enhancement, <jd oc {^ t T J t ^ , 
which physically corresponds to an Aslamazov-Larkin 
(AL) paraconductivity effect (T^]. Notice that an is more 
singular than the AL conductivity oal oc t ^ c t . This 



situation is similar to the Hall effect in superconductors 
17, 18 since both the drag and Hall conductivities are 
sensitive to the electron-hole asymmetry and require an 
expansion of L in Tji, while the AL conductivity is not. 
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Further away from the transition, the vertex suppression 
due to the pseudogap, where a oc A/T, f3 oc (A/T) 2 and 
f(T) ~ 1, results in a moderate decay of od oc , 2 v . 

Screening effects on the Coulomb potential provide an 
additional suppression of <jd for large separation between 
the layers. 

In the case when both layers are superconductors, we 
use two vertices from Eq. (161 in Eq. Q and find 
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where X ss (d) = ( ^$(fj d) > , ffCO = {^f 
with shorthand notation {a, /3}i,2 = ft}(js/ ^1,2)- 
This result shows a similar but potentially more singu- 



lar behavior near the transition as Eq. (17), however, 



with a much faster power-law decay at higher temper- 



atures, 0~ Y) oc 



since g(T) oc (AiA 2 /T 2 ) for 



T 2 ln 5 T] 

T > max{T c i, T C 2}. This high-T suppression mechanism 
due to the pseudogap was recently proposed as an expla- 
nation of the fast fall-off of the Nernst signal in under- 
doped cuprates [ID] . 

Summary- In conclusion, we developed a theory of a 
Giaever transformer based on bi-layers of cuprate super- 
conductors, assuming that the leading source of recti- 
fication originates from fluctuating Cooper pairs. The 
coupling between the layers is provided by a screened 
Coulomb interaction. The drag signal shows a strong 
singularity near the superconducting transition due to 
the paraconductivity effect, while the pseudogap leads to 
a faster drop-off at higher temperatures. Such a char- 
acteristic temperature dependence may distinguish drag 
due to preformed pairs from that due to vortices. In par- 
ticular, Ref. |2"U1 predicted that the drag resistivity po in a 
normal-superconductor bilayer due to the inductive cou- 
pling effect of vortices should follow the flux flow conduc- 
tivity <7„ of the superconductor, pn oc (h/e 2 ) 2 a v . Clearly, 
this result is very different from the contribution to pp 



originating from fluctuating Cooper pairs [Eq. (17)]. The 
general condition to have a vortex drag effect is to have 
an inhomogeneous magnetic field in the barrier region 
(otherwise, the vortices between the two layers cannot 
'find' each other). This means that the vortices must be 
well defined. As a consequence, any vortex drag regime 
above T c should have a narrow temperature range. This 
is especially true for cuprates, given the weak coupling 
of vortices between layers (i.e., pancake vortices). That 
coupling gets even weaker for undcrdoped compounds, 
which is exactly the regime we are interested in. 
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